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Abstract. We investigate the role of dust in star formation activity of extremely metal-poor blue compact dwarf 
galaxies (BCDs). Observations suggest that star formation in BCDs occurs in two different regimes: "active" and 
"passive". The "active" BCDs host super star clusters (SSCs), and are characterised by compact size, rich H2 
content, large dust optical depth, and high dust temperature; the "passive" BCDs are more diffuse with cooler 
dust, and lack SSCs and large amounts of H2. By treating physical processes concerning formation of stars and 
dust, we are able to simultaneously reproduce all the above properties of both modes of star formation (active 
and passive) . We find that the difference between the two regimes can be understood through the variation of the 
"compactness" of the star-forming region: an "active" mode emerges if the region is compact (with radius S~ 50 
pc) and dense (with gas number density J> 500 cm -3 ). The dust, supplied from Type II supernovae in a compact 
star-forming region, effectively reprocesses the heating photons into the infrared and induces a rapid H2 formation 
over a period of several Myr. This explains the high infrared luminosity, high dust temperature, and large H2 
content of active BCDs. Moreover, the gas in "active" galaxies cools (^ 300 K) on a few dynamical timescales, 
5_l ■ producing a "run-away" star formation episode because of the favourable (cool) conditions. The mild extinction 

and relatively low molecular content of passive BCDs can also be explained by the same model if we assume a 
diffuse region (with radius J> 100 pc and gas number density 5s 100 cm -3 ). We finally discuss primordial star 
formation in high-redshift galaxies in the context of the "active" and "passive" star formation scenario. 
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1. Introduction Recently Hirashita & Ferrara l|2(J(J21 hereafter HF02) 

. have proposed that the dust enrichment in extremely 

The surfaces of interstellar dust grams are known to be , , . , , . , . . 1 c ,1 1 

metal-poor primeval objects is essential lor the enhance- 

sites where efficient formation of hydrogen molecules (H2) , <■ , c ,■ ■, ,m ■ • , 

■ tt r i-i merit of star formation activity, Iheir scenario suggests 

takes place. Without dust grains, H9 forms only m the , , f , , , , , , j^r 

, . _ 1 , the following cycle between dust production and star tor- 
gas phase with a production rate much lower than the , . /M . , , . ,. „ TT 
° . . , _ , , . _. , , mation: (1) massive stars end their lives as lype 11 super- 
dust surface reaction (e.g., Peebles & Dickc 1968; Matsuda /Q AT TT\ 1-1, 111 • J. 1 

v ■ ' ; novae (SJNe 11), which supply dust grains into the mterstel- 

et al. The shielding of H 2 dissociating photons by ^ medium (IgM) (Kozaga ^ ^ ^ Todini & Ferrar& 

dust grains also enhances mo ecule formation. In general ^ Nozaw& et ^ ^ Schneider et aL[2 )04); (ii) those 

dust also absorbs a part of the radiation from stars and . , r ,• r , , 1 j • u- 1, 

f , grains enhance the formation of molecular clouds m which 

reemits it m the infrared (IR) , thereby modifying the , c c , , , . , 

. y ,' ' . r ' , . , stars lorm; (in) some ol those stars are massive and the 

spectral energy distribution SFD ol galaxies (e.g., Suva , . , , , . . r™. 1 /1 

rrrr7mi i 11 1 °- ust su PPby described m (1) occurs again. I his cycle (ij- 

et al. .1998). Therefore, dust plays an important role m a ,. , , , r ,. , , 

~ * . . , . ( mi) effectively enhances the star formation rate as shown 

both chemical and radiative properties of galaxies. ^ HF02 
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„ -i ... _ , . In fact, a large amount of dust has been suggested to 

c-mau: nirasnitaOu.pnys .nagoya-u. ac . jp ' 00 ,„,„. 

* Research Fellow of the Japan Society for the Promotion of exist at hi g h redshift (high z) (e.g., Small et al. HES). 

Science. However, it is not easy to explore the first dust enrich- 

1 In this paper, IR indicates the wavelength range where the ment in primeval galaxies at high z (z £ 5) with the 

emission from dust dominates the radiative energy (roughly present observational facilities. Therefore, nearby tem- 

8-1000 fj,m). plates for primeval galaxies are useful to test galaxy for- 
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mation scenarios. The best candidates for such a template 
are metal-poor blue compact dwarf galaxies (BCDs), since 
they are at the initial stage of metal enrichment and their 
current star formation activity is generally young (Searle 
& Sargent H1Z21 Kunth & Ostlin 127)071 . In other words, 
BCDs can be used as laboratories in which to study high-z 
primeval galaxies. 

Two classes of star formation activity have recently 
emerged observationally, as proposed for a BCD sample by 
Hunt et al. p)03al hereafter HHTIV: l27)04b)> . They argue 
that the star-formation activity in the two most metal- 
poor galaxies, SBS 0335— 052 and IZwl8, shows very dif- 
ferent properties, in spite of their similar metallicities 
(1/41 and 1/50 Zq, respectively; Skillman & Kennicutt 
1993] Izotov et al. I1999|l . The major star-forming region 
of SBS 0335—052 is compact and dense (radius r$F ^ 40 
pc, number density n J> 600 cm -3 ; Dale et al. 120011 Izotov 
& ThuanEnH. Moreover, SBS 0335-052 hosts several su- 
per star clusters (SSCs), detectable H2 emission lines in 
the near- infrared (Vanzi et al. 12000(1 . a large dust extinc- 
tion (A v ~ 16 mag; Thuan et al. lTM?! Hunt et al. lMFTl 
Plante & Sauvage 2002 ) , and high dust temperature (Hunt 
et al.EOnil Dale et al. l2OTTl Takeuchi et al. E?UU3|l . On the 
contrary, the star-forming regions in I Zw 18 are diffuse 
(tsf <i 100 pc, n Si 100 cm -3 ), and contain no SSCs. Near- 
infrared H2 emission has not been detected (Hunt et al., 
private communication), and the dust extinction is mod- 
erate (Ay ~ 0.2 mag; Cannon et al. I2002f) . We call a re- 
gion with such properties "passive" following HHTIV. The 
similar metallicities of SBS 0335— 052 (active) and IZwl8 
(passive) imply that the chemical abundance is not a pri- 
mary factor in determining the star- forming properties. 
We argue that the compactness of star- forming regions, 
which affect gas density, gas dynamics, and so on, is im- 
portant in the dichotomy of active and passive modes. 

Hirashita et al. EEEI (hereafter HHF02) show that the 
IR luminosity, the dust mass, and the rich H2 content 
of SBS 0335—052 can be explained through dust accu- 
mulation by successive SNe II in a compact (tsf Si 100 
pc) region. Moreover, SBS 0335—052 is not unique among 
BCDs; the BCDs with dense compact star-forming regions 
similar to SBS 0335—052 were dubbed "active" , and tend 
to be characterized by high surface brightness (see Figure 
1 of HHTIV). The physical state of ISM is also similar 
among "active" BCDs. They always have large dust ex- 
tinction, and a significant fraction of stellar radiation is 
reprocessed into IR; the dust temperature is also high 
and there is a "hot" component with 600-1000 K (Hunt 
et al. I2002f) . Dust properties can be further constrained 
by future IR observations for a large sample of BCDs 
(e.g., Takeuchi et al. l2T?0ljl . HHTIV also demonstrate that 
there are BCDs with converse properties, namely "pas- 
sive" ones, which are more diffuse, less dense, and of lower 
surface brightness. A representative of this category is 
I Zw 18. Contrary to "active" BCDs, IR dust emission has 
not been detected so far in "passive" BCDs. Although 
their star formation rate is lower than "active" BCDs, 
"passive" BCDs are actually forming stars, and are com- 



pletely different from passively evolving galaxies such as 
ellipticals. 

The above clear difference in dust properties implies 
that in addition to the compactness, dust should be con- 
sidered as a key to understand the "active" and "passive" 
modes, and hence to understanding how star formation 
occurs in extremely low-metallicity environments. We con- 
sider the role of dust in various compactness of regions by 
using our theoretical framework. This paper is organised 
as follows. First, in Section we explain the model that 
describes the evolution of dust content and gas state in 
a star- forming region. Then, in Section [21 we present our 
results concerning the differences between "active" and 
"passive" star-forming regions. In Section we discuss 
the interpretation of observational properties of active and 
passive BCDs in the context of our model, and consider in 
particular the two prototypes, I Zw 18 and SBS 0335—052. 
Implications for high-redshift star formation are described 
in Section and our conclusions are presented in Section 

m 

2. Model description 

For our calculation of dust amount in a star-forming re- 
gion, we essentially use the model by HF02. We treat the 
chemical reaction network concerning H2 formation in a 
way consistent with dust amount evolution. We approx- 
imate star-forming regions as homogeneous spheres and 
adopt a representative value for each physical quantity 
(i.e., we treat a star-forming region as one zone). 

2.1. ISM evolution 

The star formation process is affected by the physical state 
of the ISM. In metal-poor environments, cooling by molec- 
ular hydrogen plays an important role in star formation 
(Galli & Palla lTM8l Abel et al. 12770711 Bromm et al. 127)011 
Nishi et al. ITM81 Omukai 127)07)1 Kamaya & Silk 127)1121 
Ripamonti et al. I2002J) . The abundance of H2 should be 
ultimately a key parameter since stars are only seen to 
form in molecular complexes (see e.g., Wilson et al. [S000 
for a recent observation of a nearby galaxy). Moreover, 
enhancement of molecular gas formation is shown to re- 
sult in an active star formation (for recent results, see 
e.g., Walter et al. 2002J), and indeed there is a correlation 
between molecular amount and star formation rate (e.g., 
Bendo et al. 2002). We define the molecular fraction, /h 2 , 
as 

/h 2 = 2n H2 /n H , (1) 

where «h and nn 2 are the number densities of hydrogen 
nuclei and hydrogen molecules, respectively (i.e., if all the 
hydrogen nuclei are in the molecular form, /h 2 = 1). 

We then calculate the time evolution of ionisation de- 
gree (x), molecular fraction (/h 2 )> and gas temperature 
(T) of the hydrogen gas. The helium content is neglected 
here, since its effect should be negligible for the three 
quantities (Kitayama et al. 2001, HF02). Therefore, the 
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number density of the gas (n) is approximated by the hy- 
drogen number density (i.e., n ~ ^h)- The ionisation de- 
gree and the temperature affect the formation rate of H2 . 
The processes are treated as follows (see HF02 for details) . 

We calculate the time evolution of the ionisation de- 
gree x by taking into account collisional ionisation, re- 
combination, and photo-ionisation of hydrogen. For the 
photo-ionisation, we have included the radiative transfer 
effect following the simple recipe in Appendix of Kitayama 
& Ikeuchi (2005, w ^o derive an analytical expression un- 
der the assumption of a power law spectrum of incident 
photons. For the spectral shape parameter a of the radi- 
ation field, we adopt a = 3, because the resulting ratio 
between H2 destroying photons and H ionizing photons is 
similar to realistic OB stars. Kitayama et al. I|2001[) have 
examined a = 1 and a = 5. Their table 1 indicates that 
a = 1 overproduces the ionizing photons while a = 5 over- 
produces the H2 destroying photons. Therefore, we adopt 
an intermediate value: a = 3. In any case, different values 
of a do not affect the behaviour of x, /h 2 , and T as func- 
tions of time, because the dust accumulation (which does 
not depend on a) is the greatest influence on these three 
quantities. The normalisation of the intensity, /o(^hi), is 
determined from 



L 



cuuv 



uv 



(r SF /2) 2 



= 47T 



h{vm) ( — ) dv 



(2) 



where i/ m ; n (~ 10 15 Hz) is the minimum frequency where 
OB stars dominate the radiative energy, wtjv is the "typ- 
ical" interstellar radiation field estimated at the half of 
the radius (rgp/2; half is for a rough average), z^hi is the 
ionisation frequency of neutral hydrogen (vm = 3.3 x 10 15 
Hz) (see HF02 for details), and c is the light speed. 

In calculating /h 2 , we consider H2 formation both in 
the gas phase (via H~ or H^, where the latter is negligi- 
ble), and on the dust grain surface; H 2 destruction occurs 
through collisions with H + , H, and e~, and photodissocia- 
tion. For self-shielding effects which prohibit H2 dissocia- 
tion, equation (17) of HF02 is used, but here we substitute 
raise with tsf/2. For the reaction rates, we adopt table 1 
of HF02, but for the reaction on dust grains, i?dust, we use 
the following expression which includes the dependence on 
a: 



' 2.8x10-- (io^k) 172 ^ 



R 



a 

03 /im 



dust 



2 g cm 



cm 3 s- 1 if T < 300 K (3) 



I if T > 300 K 



In order to calculate the temperature evolution, cool- 
ing and heating must be considered in our model. We as- 
sume that cooling processes comprise collisional excitation 
and ionisation of atomic hydrogen (when T <; 8000 K) 
and collisional excitation of molecular hydrogen (which 
dominates for T < 8000 K). For the heating by stellar 
UV radiation, we have included the radiative transfer ef- 



fect by following a simple analytical recipe in Appendix 
of Kitayama & Ikeuchi ( 2000 ) . 

2.2. Star formation rate 

Stars form as a result of the gravitational collapse of a gas 
cloud. Therefore, it is physically reasonable to relate the 
star formation rate with the free-fall timescale of gas. We 
consider a star-forming region with a gas number density 
of n ~ njj- The free-fall time, tg, is estimated as 



1 x 10 7 



100 cm" 3 



-1/2 



yr, 



(4) 



where p — mn?iH (tih is the mass of a hydrogen atom) is 
the mass density of the gas. The star formation rate, is 
estimated as 



gas 



Ul 



9(t) 



0.1 



gas 



V 0.1 / V10 7 M 

M Q yr" 1 , 



(w^r e «> 



(5) 



where M gas is the total gas mass (both molecular and 
atomic) of the star-forming region, esf is the star forma- 
tion efficiency defined as the conversion efficiency of a gas 
into stars over a free-fall time, and <d(t) is Heaviside's 
step function [Q(t) = 1 if t > and 6(t) = if t < 0]. 
Thus, we define the zero point of time t at the onset of 
star formation in the star-forming region. For simplicity, 
we assume a constant star formation rate (i.e., csf, M sas , 
and tih are approximated to be constant) and a spheri- 
cal homogeneous star forming region. In reality, there is 
probably a significant density inhomogeneity; however, the 
star formation rate in the entire star-forming region may 
be regulated by the dynamical time governed by a spa- 
tially averaged density, and the homogeneous density can 
be regarded as such an average. 

The hydrogen number density can be related to gas 
mass as 



4"7T 

— r| F n H m H = Mgas , 



(6) 



where tsf is the radius of the star-forming region. Thus, 
the numerical value of the number density is estimated as 



n H ^ 100 



^SF 

100 pc 



Af„. 



10 7 M Q 



We also define the gas consumption timescale t gas : 
M„ 



t. 



gas tff 



(7) 



(8) 



Since the effect of gas conversion into stars becomes sig- 
nificant for t ~ i gas , we stop the calculation at t = i gas /2. 
After this time, star formation may be suppressed because 
of the lack of gas. 
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2.3. Evolution of dust content 

We calculate the evolution of dust mass in a star-forming 
region considering the dust supply from SNe II. The pre- 
existing dust at t — in the star-forming region is ne- 
glected here. This is a good approximation for the ex- 
tremely metal-poor BCDs treated in this paper. In fact, 
a large fraction of dust amount contained in some BCDs 
can be explained by the present star formation activity 
(see Section 13. 3f) . Preexisting dust would strengthen all 
the dust effects described below. 

If the starburst is young, and there has been no pre- 
vious burst, dust can be supplied only from massive stars 
with short lifetimes. Thus, SNe II are the dominant source 
of dust formation if the age is typically less than a few 
times 10 s yr (Dwek [TSHSJ) . Thus, the dust formation rate 
is related to the rate of SNe II. The SN II rate as a function 
of time, 7(t), is given by 



lit) 



i (f>(m) dm , 



(9) 



8 A/(7 



where ij)(t) is the star formation rate at t (Eq.JSJ, <fi(m) is 
the initial mass function (IMF; the definition of the IMF 
is the same as that in Tinslev liySOII . r m is the lifetime of 
a star whose mass is m, and we assume that only stars 
with m > 8 Mq produce SNe II. In this paper, we assume 
a Salpeter IMF {4>(m) cx m -2 ' 35 ) with the stellar mass 
range of 0.1-60 Mq (HF02). If we assume a higher upper 
limit such as 100 Mq, we expect a higher SN II rate and 
dust production rate. However, if the progenitor mass m 
is larger than ~ 50 M©, almost all the produced metals 
may collapse into a central blackhole (Tsujimoto et al. 
I1997f> . Todini & Ferrara 1)200 1[1 treated only progenitor 
stellar mass less than 40 Mq, and the dust production 
rate for massive stars (<; 60 Mq) is unknown theoretically. 
Recently, however, there have been some advances in this 
field (Nozawa et al. l2o"0"51 Schneider et al. l2"0"0l|l . and we 
are developing a model to include the effect of massive 
stars. 

Dust destruction by shocks from SNe II can be im- 
portant (Dwek & Scalo HMH Jones et aL fiMfy . The de- 
struction timescale tsn is estimated to be (McKee 1989 
Lisenfeld & Ferrara IT5M)> 



M 



tsn 



gas 



(10) 



where M s is the mass accelerated to a velocity large 
enough for dust destruction by a SN blast (~ 100 km s _1 ), 
7 is the SN II rate, £sn is the efficiency of dust de- 
struction in a medium shocked by a SN II. We adopt 
M s = 6.8 x 10 3 Mq (Lisenfeld & Ferrara [15551 see also 
Tielens HH3> and e SN = 0.1 fMcKee ll"5g5|) . Since we are 
interested in extremely metal-poor environments, we as- 
sume the relation between stellar mass and lifetime of 
zero-metallicity stars in table 6 of Schaerer (2002). 2 



The case without mass loss is applied for consistency with 
HF02. This does not unduly affect our results however since 



Dust grains can also be destroyed by strong and hard 
interstellar radiation fields (e.g., Boulanger et al. [1988 
Puget & Leger 115551 Voit 1T5521 Contursi et al. JMty . 
However, the grain populations most strongly affected are 
the carriers of the aromatic band features (e.g., Polycyclic 
Aromatic Hydrocarbons) and very small grains (diame- 
ter a < 0.01/^m), neither of which are considered here. 
Hence, since grain destruction processes by UV radiation 
are most severe for very small grains, and in any case are 
poorly known, we will neglect them. 

In this case the rate of increase of Md us t is written as 



dM dust M dust 
■ «i dust 7 - 



dt 



(11) 



where mdust is the typical dust mass produced in a SN II. 
Todini & Ferrara (|2001|1 show that mdust varies with pro- 
genitor mass, metallicity (Z), and input energy of SNe II. 
The Salpeter IMF-weighted mean of dust mass produced 
per SN II for (1) Z = 0, Case A, (2) Z = 0, Case B, (3) 
Z = IQ- 2 Zq, Case A, and (4) Z = 10~ 2 Z Q , Case B are 
(1) 0.22 Mq, (2) 0.46 M Q , (3) 0.45 M Q , and (4) 0.63 Mq, 
respectively [Z is the metallicity, and Cases A and B cor- 
respond to low and high explosion energy 3 for Case A and 
Case B, respectively.) We adopt the average of the four 
cases, i.e., mdust = 0.4 Mq (HHF02). This may not be an 
overly optimistic estimate given the low upper mass cut- 
off of 35 Mq adopted by Todini & Ferrara i|2770T|l . and the 
higher dust production estimates of Nozawa et al. (2003). 
The calculated dust mass Md us t is roughly proportional 
to mdust) and the timescale on which the dust extinction 
effects appear is approximately proportional to l/m d ust- 

2.4. Evolution of metal content 

The evolution of metal content can be calculated once 
the star formation history and metal yield are fixed (e.g., 
Tinslev H980fl . Because of the young age range (^ 10 8 yr) 
treated in this paper, we assume that the metal produc- 
tion is dominated by SNe II. The evolution of the mass of 
metals (a given species is indicated by i) in the gas phase, 
Mi, is calculated by subtracting the dust mass from the 
metal mass. This is expressed as 



dM t 
dt 



dM, 



ma - 



dust, 



dt 



(12) 



where m^ is the IMF-weighted average of metal mass (for 
species i) ejected per SN II and M du st, t 1S the mass of 
element i in dust phase. We calculate the abundances of 
carbon and oxygen, because those two elements are the 
two major metals produced in SNe II (Woosley & Weaver 
IT555T) . Following HF02, we adopt M dus t,o = 0.15M du st, 
M d ust, c = 0.36M du st, m = 1.2 Mq, and m c - 0.17 M Q ; 
i.e., we assume that carbonaceous grains are responsible 



the differences in metallicity and luminosity for the case with 
and without mass loss are within a factor of 2. 

3 The kinetic energies given to the ejecta are 1.2 x 10 51 erg 
and 2 x 10 51 erg 
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for 36% of the total dust mass, and the other grains con- 
tain oxygen with mass fraction of 23% (i.e., the oxygen 
fraction in mass is 0.64 x 0.23 = 0.15). 

2.5. Radiative properties 

One of the most direct ways to constrain dust content is to 
observe IR continuum emission from dust grains. We now 
derive the evolution of IR luminosity. Because of the large 
cross section of grains against UV light and the intense UV 
interstellar radiation field in star-forming galaxies, we can 
assume that the IR luminosity is equal to the absorbed 
energy of UV light. 

2.5.1. UV and IR 

First, we estimate the fraction of the UV radiation ab- 
sorbed by dust. We define the following typical optical 
depth, T dust , as 



Tdust = 7ra 2 Quv"-dust?'SF/2 . 



(13) 



where 7ra 2 QTjv is the typical absorption cross-section for 
UV light (Quv is the dimensionless absorption cross- 
section normalised by the geometrical cross-section), a is 
the dust radius, ridust is the mean dust grain number den- 
sity, and we assume that the typical path length for a UV 
photon is half of rgp as a rough average. 4 In this paper 
we assume for simplicity a single value for a, consistently 
with the typical size of dust grains produced by SNe II 
(e.g., Todini & Ferrara I2TM1 Nozawa et al. EDD^- The 
dust mass density 6 is related to the mean dust number 
density n dust as 



in 

T 



in 



a 3 Sn dust —r? 



SF 



= Af dl 



where S is the grain material density. By solving Eq. (|14fl 
for ndust and substituting it into Eq. H13[l. we obtain 



Tdust 



9 QuvM, 



dust 



32tt 



iSr, 



(15) 



SF 



We define the "attenuation function" -E(rd us t) as the 
fraction of UV light escaping from the star-forming region. 
Thus, we write the UV luminosity of the star-forming re- 
gion as 



^uv 



(16) 



where Ltjv.o is the intrinsic UV luminosity. The functional 
form of -E(rdust) depends on the geometry of dust distri- 
bution, and we examine the following two representative 
simple cases. One is the "screen" distribution, in which all 
the dust is located at the radius of tsf from the centre of 



4 In order to fully derive the precise path length instead of 
roughly dividing tsf by 2, a detailed modelling of radiative 
transfer including dust scattering is necessary. Here we only 
mention that with a given dust mass, Td us t linearly scales with 
the adopted path length. 



the star-forming region. The other is the "mixed" geom- 
etry, in which the spatial distributions of dust and stars 
are the same (i.e., the mass ratio between stars and dust 
is constant). In the screen geometry, 



E(r dv 



E sc (Tdu 



exp(- Tdust) , 



(17) 



while in the mixed geometry, 



E(T dust ) = E mix (T dust ) = [1 - exp(-T du st)]/Tdust ■ (18) 

Since E ml ^ > E sc for a given Tdust, the screen geometry 
shields the UV light more efficiently. Moreover, the ex- 
ponential behaviour of E sc results in a strong cut-off of 
UV light at a certain time when a significant amount of 
dust is accumulated, while E mlK ~ 1/Td us t for Td us t 3> 1, 
which means that UV light originating from the "surface" 
always escapes. We do not treat a clumpy dust distribu- 
tion, since inhomogeneity changes average optical depths 
in ways that depend on the specific geometry (Natta & 
Panagia 1984J; such a treatment is beyond the scope of 
this paper. 

-^uv,o is assumed to be equal to the total luminosity 
of OB stars, whose mass is larger than 3 M Q (Cox [2000): 



^uv,o(*) 



dm 



3 Ma 



dt' L(m) 4>{m) t/j(t - t') , (19) 



where L(m) is the stellar luminosity as a function of stellar 
mass (to). This UV luminosity is also used in Eq. (|16|l to 
estimate the UV luminosity Ltjv after the dust absorption, 
and the same Ltjv is also in Eq. (J2J. For £(m), we adopt 
the model of zero-metallicity stars without mass loss in 
Schaerer ( 2002). Assuming that all the absorbed energy is 
rccmittcd in the IR, the IR luminosity Lir becomes: 



(14) L m = L 



UV,0 



£tjv — L 



UV,0 



1 - E(T dust )} 



(20) 



2.5.2. Dust temperature 



Another important quantity representative for the radia- 
tive properties of a star-forming region is dust temper- 
ature. The equilibrium dust temperature is determined 
from the balance between incident radiative energy and ra- 
diative cooling. The equilibrium temperature is expressed 
as follows (Takeuchi et al. 



T ~- ( hc \ { 63m uvQuv \ 1/6 
dust ~ \nk) \ Mir{2irAa)hc J ' 



(21) 



where the UV radiation field is calculated by Eq. (J2J) and 
the dimensionless dust emissivity in the IR is assumed to 
satisfy 5 



(a, A) 



2nAa 



(22) 



For silicate grains A = 1.34 x 10 3 cm (Drapatz & Michel 
I1977|) . while for carbonaceous grains A = 3.20 x 10 -3 cm 



5 Strictly speaking, this is valid for A <; 20 fj,m for silicate 
grains, and A <; 50 fim for carbonaceous ones. 
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(Draine & Lee ITQgl Takeuchi et al. 12005|). The deriva- 
tion of Eq. (|21() is found in the Appendix. The numerical 
expression for Td us t is 



l dust 



14.6 



mtjv 



4.0 x 10~ 14 era: cm 



1/6 



.4 



3.2 cm x 10" 



-1/6 



0.03 /zm 



-1/6 



where we explicitly express the dependence on a. The equi- 
librium dust temperature is seen to be rather insensitive to 
changes in the input parameters. We adopt A = 3.2xl0~ 3 
cm since that is the grain radius appropriate for graphites 
(Todini & Ferrara 2001). With this value, we obtain 



'IB 



(a, A) = 6.03 x 10" 



0.03 J V!00 A*m 



A 



(24) 



For a grain size of 0.03 Atm (carbonaceous grains), and 
a density 6 = 2gcm~ 3 (Draine & Lee I1984*)l . this cor- 
responds to the following value of dust mass-absorption 
coefficient, n v (= 3Qm/4aS) (Hildebrand 19S5J): 



75 



A 



100 /im 



2 g cm 



2 —1 

cm g 



(25) 



Although there is little observational constraint on Qir 
and k v of dust produced by SNe II, we can compare our n v 
with the data of Galactic or extragalactic observations. At 
A = 850 A*m, k u ~ 1.0 cm 2 g _1 , while observational data 
suggest k v ~ 0.7-2.4 cm 2 g" 1 (Alton et al. 120011 James 
et al. 12002)) . Bianchi et al. l(T000)) find that Q m around 
A = 100 Aim is roughly ~ 10 _3 -5 x 10~ 4 if we assume Q 
at V band is ~ 1 (see also Fig. 5 of James et al. 120021 and 
references therein). 

2.6. Initial conditions 

For the onset of star formation, the gas must be cool. The 
molecular hydrogen cooling eventually cools gas down to 
T ~ 300 K (e.g., Tegmark et al. 1997). During the phase of 
efficient molecular formation in gas, the ionisation fraction 
x is roughly ~ 10 -4 and the molecular fraction is /h 2 ~ 
10~ 3 . After this phase, the ionisation degree x decreases 
and the gas phase reaction stops. Since stars begin to form 
during molecular formation, we set the initial values as 
x = 10- 4 , T = 300 K, and /h 2 = 10~ 3 . However, the gas 
starts to be ionised and the molecules begins to dissociate 
soon (< 1 Myr) after the onset of star formation, and the 
results do not depend strongly on the initial conditions for 
x and /h 2 ■ The initial temperature does not change the 
evolution of gas as long as T < 10 3 K initially. 



2.7. Selection of parameters 

For the UV light, because of the short wavelength, we can 
assume that Qtjv — !• F° r the grain radius, we examine 



the following two cases: small dust (a = 0.03 /im) proposed 
by Todini & Ferrara (2001), and large dust (a = 0.1 Aim) 
possibly produced even in SNe II (Nozawa et al. 2003). 
We fix m du st = 0.4 Mq (Section [^J, and adopt S ~ 2 g 
cm~ 3 (Draine & Lee 1984} . As mentioned in Section I2"51 
the timescale of dust effects (enhancement of molecular 
content, enhancement of UV shielding, etc.) is roughly 
inversely proportional to mdust- 

The free parameters which remain to be determined 
(23) are rsF, M sas , and egp. Since our aim is to investigate the 



J gas j 

properties of "active" and "passive" star-forming regions, 
which correspond to compact and diffuse regions, respec- 
tively, we should examine various tsf- In order to concen- 
trate on the effect of compactness, the gas mass is fixed. A 
representative value for a star-forming region in BCDs can 
be taken as M gas = 10 7 Mq. This is a typical value for gas 
masses derived for star-forming regions in SBS 0335—052 
and IZwl8 representative extremely metal-poor BCDs 
( Section 14. 2|) . In order to concentrate on the variation of 
tsf, we therefore fix M gas = 10 7 Mq, and examine a re- 
gion size range typical of BCDs: rgp = 30 pc, 100 pc, and 
300 pc, corresponding to the density n ~ 3000 cm~ 3 , 100 
cm~ 3 , and 3 cm~ 3 , respectively, with M gas = 10 7 Mq. 
The most compact case is representative of the "active" 
class, while the most diffuse case is typical of the "passive" 
class. 

The remaining free parameter is the star formation ef- 
ficiency £sf 7 which we assign a value 6sf = 0.1 (e.g., Ciardi 
et al. lMTUl Ferrara et al. 1217001 Inoue et al. l2lTOI Barkana 
I2UU21 Scannapieco et al. 1211351 Salvaterra & Ferrara 12005) 1 
unless otherwise stated. Roughly speaking, the timescale 
on which the dust effects appear scales inversely with egp . 
For example, the rapid increase of molecules and drop of 
temperature appear on a timescale proportional to 1/esF- 
It is possible that e§F may change significantly as a func- 
tion of time, but our simple assumption of constant esf 
is a useful first-order approximation, csf and the star for- 
mation rate should be regarded as an average over the 
timescale in which we are interested. High and low star for- 
mation efficiencies are investigated in the models of I Zw 18 
and SBS 0335-052 in Section IO 

These values are all representative of metal-poor star- 
forming regions. Typically, BCDs have ip ~ 0.01-1 
M Q yr" 1 , r - 0.1-3 kpc and n ~ 10-1000 cm" 3 (Popescu 
et al. lTOOOl Hopkins et al. l2002)l . 

3. Results 

Here we illustrate the results for the time evolution of 
star-forming regions. First we adopt a = 0.03 /im and the 
screen dust distribution. Then, in Section 13.41 we inves- 
tigate the dependence on grain size by assuming larger 
grains (a = 0.1 /im) and on geometry by examining the 
effects of a mixed dust distribution fEa.llSfl. 

For continuous star-forming activity, the gas should 
continue to collapse. With a fixed density, the typical mass 
for the collapse, i.e., the Jeans mass, is determined by the 
gas temperature. Therefore, the evolution of the gas tern- 
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perature is one of the most important processes for under- 
standing the fate of a star-forming region. In the following, 
we show the time evolution of gas temperature (T) and 
molecular fraction (/h 2 ) f° r various region radii. Because 
we have assumed fixed values for M gas , the density «h is 
thus determined by setting tsf- Moreover, since we have 
also fixed csf, the resulting star-formation rate ip, given 
by Eq. JSJ, also depends on rg-p through the density n 
dependence. 

3.1. Gas state 

In Fig.QJ we show the evolution of (a) gas temperature, (b) 
ionisation degree, and (c) molecular fraction for tsf = 30, 
100, 300 pc (solid, dotted, and dashed lines, respectively). 
We stop our calculation at t = i gas /2, when our assump- 
tion that M gas is constant becomes invalid (see Eq. |HJ). 
Therefore the lines are truncated at 8 Myr and 50 Myr for 
rsp = 30 pc and 100 pc, respectively. We observe that 
as the region size becomes smaller (i) the gas remains 
cooler, (ii) the ionisation degree remains lower, and (iii) 
the molecular fraction increases more rapidly. 

In the time evolution of these three parameters, accu- 
mulation of dust grains plays a fundamental role. The dust 
accumulation in a compact region results in a large optical 
depth and a large density of dust grains. As a result, the 
molecular formation tends to increase rapidly, because of 
the large density of dust grains and the consequent large 
shielding of dissociating photons. Since the UV photons 
are efficiently shielded by dust grains, the gas is kept cool 
and the ionisation degree remains low. 

The temporal behaviour of the molecular fraction is ex- 
plained as follows (see Fig.^). The first drastic increase 
of the molecular fraction, seen around 4 and 13 Myr for 
?"SF = 30 and 100 pc, respectively, is due to the activation 
of the gas-phase reactions caused by the increase of ion- 
isation degree and the temperature drop. The molecular 
formation is suppressed for tsf = 300 pc because the gas 
temperature remains high ~ 10 4 K. The second drastic 
increase of the molecular fraction, seen around 5 and 20 
Myr for tsf = 30 and 100 pc, respectively (for rsF = 300 
pc, the second increase is not before 100 Myr), is due to 
the onset of molecule formation on dust grains: dust grains 
shield the ionising photons and facilitate the gas cooling 
down to 300 K. At this temperature, the H2 formation on 
dust grains becomes possible. For rgF = 300 pc, represen- 
tative of the passive mode, H2 formation occurs only in 
the gas phase with a low reaction rate because of the high 
gas temperature. 

3.2. Luminosities 

To see how much UV light is shielded by dust, we show 
the evolution of IR and UV luminosities in Fig. for 
rsp = 30 pc and r$-p = 300 pc, corresponding to the ac- 
tive and passive modes, respectively. We see that the IR 
vs. UV luminosity ratio increases more rapidly in the com- 



pact (rgF = 30 pc) region than in the diffuse (rgp = 300 
pc) one. This is due to the large dust optical depth in 
compact star-forming regions. The efficient shielding of 
UV in compact regions also keeps the molecular fraction 
high since UV dissociation is suppressed. 

The evolution of dust temperature determined from 
Eq. (|23J) is also shown in Fig. 03 for rgF = 30 pc and 
300 pc. The lines start from the onset of the first SNe II 
(i.e., t = 3 Myr), when dust production begins. The dust 
temperature is initially ~ 100 K and decreases rapidly for 
the compact "active" case. This rapid decrease is due to 
the strong shielding of UV light by dust grains. The dust 
temperature remains ~ 50 K for the "passive" case. 

3.3. Chemical enrichment 

In order to ascertain if our chemical enrichment model is 
capable of explaining the observed metallicities, we cal- 
culate the evolution of oxygen abundance. Since the op- 
tical oxygen emission lines from ionised regions are often 
used as a tracer of chemical abundance (e.g., Izotov & 
Thuan 1999), we calculate the evolution of oxygen mass 
(see Section EU; the solar abundance is assumed to cor- 
respond to Mo /Mg as = 0.01 fCox l200T)j> . where M gas is 
assumed to be constant in time. Fig. OK shows the evolu- 
tion of [O /H] for the three region radii depicted in Fig. ^ 
The metallicity of the 30-pc region reaches 1/50 solar in 
7 Myr. Not surprisingly, the metal enrichment in diffuse 
regions proceeds more slowly: in particular, the 300-pc re- 
gion is enriched to 1/50 solar in ~ 30 Myr. 

The evolution of dust-to-gas ratio (V = Md us t/M g ) is 
shown in Fig. Et>- We have assumed that the star-forming 
regions contain no dust and metals initially. This assump- 
tion holds if a large part of the present dust and metals is 
supplied during the present episode of star formation. This 
assumption is adopted based on the picture that BCDs 
are relatively young galaxies whose major star formation 
episodes have occurred quite recently (S5 100 Myr; e.g., 
Tomita et al. 12001 Takeuchi et al. lMiljl . This point will 
be further discussed when we compare our model with the 
observational properties of BCDs fSection l4.2fl . 

3.4. Grain size and distribution geometry 

Nozawa et al. I|2003fl have recently pointed out that grains 
larger than predicted by Todini & Ferrara l|2001[l can form 
even in SNe II. Grains may be as large as a ~ 0.1 fim, and 
therefore, we also examine this case. 

In order to evaluate the effect of shielding, we show the 
evolution of UV and IR luminosities for a = 0.1 /im with 
the screen geometry in Fig.0^. With a fixed dust mass, the 
dust optical depth against the UV photons, Td us t, is pro- 
portional to a^ 1 , because ridust oc a~ 3 and Td us t oc a 2 rid U st 
(Eq. I13fl . Therefore, if the grain size is larger, more UV 
radiation escapes without being absorbed by dust grains. 
Consequently, the IR luminosity becomes dominant rela- 
tive to Ltjv later for a = 0.1 /im than for a = 0.03 /an, 
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Fig. 1. Time evolution of a) gas temperature, b) ionisation degree, and c) molecular fraction for various radii of the 
star- forming region (solid, dotted, and dashed lines for ?"sf = 30 pc, 100 pc, and 300 pc, respectively). The gas mass 
and the star formation efficiencies are assumed to be 1.0 x 10 7 Mq and 0.1, respectively. 




Fig. 2. a) Time evolution of IR and UV luminosities (solid and dotted lines, respectively). The left lines are for the 
compact (rsF = 30 pc) case and the right lines are for diffuse (r$F = 300 pc) case. The gas mass is assumed to be 
10 7 Mq (same as Fig.^l. b) Time evolution of dust temperature. The upper and lower lines are for tsf = 30 pc and 
rsF = 300 pc, respectively. The lines start from the age of around 3 Myr, when dust begins to be produced. The screen 
dust distribution and the small dust grain (a = 0.03 /im) are assumed for both panels. 



although the difference is not significant especially for the 
compact (rsF = 30 pc) region. For tsf = 30 pc, the IR lu- 
minosity dominates as soon as the dust production begins 
(~ 3-4 Myr). 

We also show the evolution of dust temperature in 
Fig. El 3 - The dust temperature drops rapidly in compact 
regions (rsF = 30 pc) because of the efficient shielding 
of UV radiation. For diffuse regions (rsF = 300 pc), the 
temperature drop relative to small grains is not significant 
because the shielding is inefficient. Owing to the a -1 / 6 de- 
pendence of dust temperature (Eq. 123(1 , the dust temper- 
ature is lower for a — 0.1 //m than a = 0.03 /im. 

Next, in order to investigate the effect of dust dis- 
tribution geometry, we examine the mixed geometry us- 
ing the attenuation function expressed in Eq. I|18fl rather 
than Eq. (|17fl . The evolution of UV and IR luminosi- 



ties is shown in Fig. [SJi. The UV light decreases more 
mildly in the mixed geometry than with the screen ge- 
ometry, because UV light originating from the "surface" 
of star-forming regions always escapes as described in 
Section 12.5.11 The evolution of dust temperature is shown 
in Fig. [3Jd. We see that the strong "exponential" drop of 
dust temperature with the screen geometry (Figs. 0d and 
0J)) is not seen in the mixed case because of the milder 
shielding of UV light. 

The physical state of the gas is also affected by the as- 
sumptions about dust grains. In order to clarify this point, 
we compare the evolution of gas temperature, ionisation 
degree, and molecular fraction in Fig. [5] for the case of 
rsF = 100 pc. The solid lines, which are the same as the 
dotted lines in Fig. ^ show the case of a — 0.03 /im and 
the screen geometry (called "standard"). The dotted and 
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Fig. 3. Time evolution of oxygen abundance in gas phase. [O/H] is the oxygen abundance in logarithmic scaling, where 
[O/H] = is the solar oxygen abundance. The solid, dotted, and dashed lines are for r$F — 30 pc, 100 pc, and 300 pc, 
respectively. The gas mass and the star formation efficiency are assumed to be 10 7 M© and 0.1, respectively (same as 
Fig.HJ. 




Fig. 4. a) Same as Fig. |2t, but for large grains (a = 0.1 /mi), b) Same as Fig. |5Jd, but for large grains. 



dashed lines represent the cases of the mixed geometry 
(with a = 0.03 /mi) and a = 0.1 /tm (with the screen 
geometry). Because the UV optical depth of dust grains 
becomes smaller as the grain size increases, the heating 
and dissociation rates are larger for larger grains, which 
explains the behaviour of the dashed lines. For the molec- 
ular fraction, the formation rate on dust grains (Eq. 
decreases because of the reduced dust surface, which also 
contributes to the slow increase of /h 2 for the larger grains. 
In the mixed geometry, the shielding of UV light is inef- 
ficient; the gas is heated and molecules are dissociated. 
Therefore, we conclude that the grain size and distribu- 
tion are important to determine the physical state of gas. 



To further investigate the effect of dust size and spa- 
tial distribution on our model results, we are planning 
to extend our study to include infrared SEDs. The study 
of SEDs has the advantage of being able to treat the 
multi-temperature, multi-size grain distributions, as well 
as stochastically heated grains (Takeuchi et al. 120031 see 
also Galliano et al. 12.003). Not only extragalactic stud- 
ies but also Galactic studies will help us understand why 
H2 formation is enhanced in actively star-forming regions 
(e.g., Habart et al. l5UU3) . 
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Fig. 5. a) Same as Fig. Et, but for the mixed geometry of dust distribution, b) Same as Fig. but for the mixed 
geometry of dust distribution. 




Fig. 6. Time evolution of a) gas temperature, b) ionisation degree, and c) molecular fraction for rg-p = 100 pc and 
Mgas = 10 7 Mq. The solid, dotted, and dashed lines represents the "standard" (a = 0.03 >tm and screen geometry) 
case, the mixed dust geometry (a = 0.03 jum), and the large dust grains of the star-forming region (solid, dotted, and 
dashed lines for tsf = 30 pc, 100 pc, and 300 pc, respectively). The gas mass and the star formation efficiencies are 
assumed to be 1.0 x 10 7 Mq and 0.1, respectively. 



4. Active and passive BCDs 

4.1. Physical properties 

In the above, we have shown that both compact and dif- 
fuse star-forming regions can keep the ISM cool and rich 
in molecules. A compact and dense (typically n 10 3 
cm" 3 , rgp ^ 50 pc, and tp <; 0.1 M & yr" 1 ) region forms 
stars on a timescale of ^ 10 Myr. This regime is termed an 
"active" mode in HHTIV. Since grains supplied by SNe II 
efficiently shield UV photons, the gas cools and molecule 
formation is enhanced. This implies that in such an ac- 
tive region, gas continues to collapse and stars form in a 
"run-away" mode. On the contrary, if a star-forming re- 
gion is diffuse (n ^ 50 cm" 3 and r$F ^ 100 pc), the star 
formation takes place "quiescently" (ip ^ 0.07 M Q yr" 1 ) 
on a timescale of > 10 8 yr. This mode is called a "passive" 
mode in HHTIV. The latter regime inhibits an "active" 



mode for the following two reasons: first, the gas collapse 
occurs slowly because of a long dynamical timescale; sec- 
ond, because of the inefficient shielding, it is subject to 
strong UV heating and photodissociation of molecules re- 
sulting from the star formation activity. 

Let us consider our results from the viewpoint of the 
gas state. In an active region, independently of the ini- 
tial conditions, the efficient dust accumulation lets the 
gas cool to T ^ 300 K through molecular cooling and UV 
shielding. Such a low temperature is crucial for the effi- 
cient formation of H2 (Cazaux & Tielens ;2002), causing 
a rapid increase of /h 2 • On the contrary, the evolution of 
a passive region proceeds without molecular cooling and 
UV shielding, and remains hotter (T 5000 K). This im- 
plies that the gravitational collapse in diffuse regions is 
inefficient and self-regulated, because of the stellar heat- 
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ing, especially the photo-ionisation heating (Lin & Murray 
1992). 

What determines the size and density of a star-forming 
region? With a constant density, the typical length of a 
gravitationally bound region is given by the Jeans length 
Aj: 



A., 



100 



T 



10 4 K 



1/2 



100 cnr 



-1/2 



pc. 



(26) 



If a star-forming region is photo-heated, the gas temper- 
ature becomes T ~ 10 4 K. Therefore, the typical size of 
a self-regulated star-forming region is determined by the 
gas density as ~ 100(n/100 cm -3 ) -1 ' 2 , consistent with 
the size of passive star-forming regions. Since an active 
star-forming region further cools through dust shielding, 
the gas collapses. We suggest that this collapse finally pro- 
duces the high surface brightness and compact size ob- 
served for active star-forming regions. 

Therefore, the initial gas density, which determines the 
self-gravitating length (Eq. I26[) of a star-forming region, 
is important for the bifurcation into active and passive 
regimes. However, ambient pressure can also influence the 
fate of a star- forming region (e.g., Elmegreen & Efremov 
(1997. Elmegreen & Elmegreen ll997l Elme green & Hunter 
2000} . Elmegreen & Efremov (|1997[1 stress that the star 
formation efficiency is affected by ambient pressure. Shock 
compression may be important to create high-pressure en- 
vironments (Elmegreen & Elmegreen 1997). A high pres- 
sure environment favours the formation of SSCs (Bekki & 
Couch USD Billett et al. EEEJ), which are observation- 
ally known to be associated with strong starbursts (e.g., 
Hunter et al. I1994J) . If a star-forming region is confined 
by a high environmental pressure, the compression of gas 
can lead to a high density and a short free-fall time. This 
means that high pressure regions tend to engender an ac- 
tive mode of star formation. Since there is no reason why 
the ambient pressure should be related to metallicity, this 
naturally explains why metallicity is not a primary factor 
of distinguishing active and passive modes. 

To link ambient pressure and region size, it is ap- 
propriate to utilise the Ebert-Bonnor formalism for self- 
gravitating, pressure-bounded isothermal spheres (Ebert 
1955 Bonnor fl956fl . The relation indicates that if rgp is 
smaller than the critical radius, 



r c = 0.49 



43 



kT ( 1 



1/2 



m H \Gp C xt, 
T \ / nT 



10 4 K / V 10 6 cm" 3 K 



-1/2 



pc, 



(27) 



the region becomes unstable and collapses. Therefore, if 
the region is compressed roughly down to 2r c ~ 100 pc, 
it becomes unstable, and can evolve into an active star- 
forming region. Another important role of ambient pres- 
sure is to avoid the expansion of the star-forming region 
and to maintain a high-density environment. In summary, 
ambient pressure has two effects: one is to make the star- 
forming region gravitationally unstable, and the other is 



to keep the density high which leads to short free-fall times 
(e.g., ElmegreenpOOO). 

4.2. Two prototypes 

One of our main motivations is to clarify the reason why 
SBS 0335— 052 and IZwl8 show different modes of star 
formation: active and passive (HHTIV). First, we "sim- 
ulate" these two galaxies with our models. The neces- 
sary quantities for our model are M gas and tsf- The Hi 
gas mass has been observationally derived as ~ 1O 9 M0 
(Pustilnik et al. IMH) and 2.6 x 10 7 M (van Zee et 
al. EMS) for SBS 0335-052 and IZwl8, respectively. 
However, we should consider these values as upper limits 
for M gas , since M gas is the gas mass in the star-forming 
regions, not in the entire system. It is difficult to resolve 
the star-forming regions in BCDs because they are gen- 
erally small. On the other hand, typical gas densities in 
star-forming regions have been derived from the argument 
of collisional excitation based on high (spatial) resolution 
spectroscopy (e.g., Izotovet al. 1999). Although this argu- 
ment is biased toward the ionised region close to massive 
stars, the gas density derived in this method can be consid- 
ered to be representative of the entire star-forming region. 
Therefore, we will attempt to distinguish active and pas- 
sive star formation and constrain the gas mass M gas on 
the basis of observed region sizes r$F (e.g., HHTIV) and 
gas densities. 

4.2.1. IZwl8 

We adopt r^F — 100 pc and n = 110 cm -3 (Izotov et al. 



1999, HHTIV). Eq. 10 indicates M gas ~ 1 x 10 7 M Q which 
is close to the M gas estimated by van Zee et al. fl998|. The 
free-fall time become tg ~ 10 Myr (Eq. 0}. The observa- 
tional star formation rate is ~ 0.04-0.1 M© yr _1 (e.g., 
Hopkins et al. 2002). This range of star formation rate 
is consistent with Eq. (JSJ) if we assume e = 0.04-0.1; we 
adopt e = 0.07 as a rough central value. In Fig. we 
show the evolution of various physical quantities typical 
of I Zw 18 by dotted lines. The shielding of UV photons 
by dust at t > 20 Myr causes the efficient cooling, recom- 
bination, and molecular formation. In particular, if the 
age of I Zw 18 is 20 Myr, the poor molecular content 
is consistent with observations (e.g., Vidal-Madjar et al. 

E3HHJ- 

Kamaya & Hirashita (2001} suggest that H2 forma- 
tion in I Zw 18 may be regulated by the gas phase reac- 
tion, rather than by the formation on dust surfaces. If 
we assume the age of the youngest burst in I Zw 18 to 
be ~ 10 Myr, our calculation gives an ionisation fraction 
of ~ 10~ 3 , roughly consistent with Kamaya & Hirashita 
(200]}. The molecular fraction is ~ 10~ 6 , consistent also 
with the upper limit by Vidal-Madjar et al. (2000). 

Recchi et al. (2002} have proposed two instantaneous 
bursts separated by a quiescent period for IZwl8. The 
age of the recent burst deduced by them (4-15 Myr) is in 
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Fig. 7. Time evolution of a) gas temperature, b) ionisation degree, and c) molecular fraction for the model of 
SBS 0335— 052 and IZwl8 (solid and dotted lines, respectively). The dashed line continuing from the solid line shows 
the evolution after i gas /2, when more than half of the initial gas mass is converted into stars. 



the range treated in our paper, but an older burst which 
occurred ~ 300 Myr is probably also necessary (see also 
Aloisi et al. lTMfl Ostlin lMTUjl . Hunt et al. (|2003b|) also de- 
rive a young burst age of 3-15 Myr, and a previous episode 
of star formation no older than < 500 Myr. However, the 
stars with the age < 10 Myr contribute to the gas heat- 
ing more than the old stars, and therefore, the age in our 
model should be taken as the age of the more recent burst. 

In Fig. |H1 we present the evolution of oxygen abun- 
dance and dust-to-gas ratio for the model of I Zw 18 (dot- 
ted lines). The horizontal dash-dotted line represents the 
metallicity level observed for IZwl8 (1/50 solar). This 
metallicity is reached at the age of ~ 20 Myr, which is 
roughly consistent with the above ages for the current star 
formation activity. At this age, the dust-to-gas ratio be- 
comes 10~ 4 . The observational constraint on the dust-to- 
gas ratio is difficult at the moment, since we must specify 
the gas mass contained in star-forming regions, which are 
generally small for BCDs, and thus difficult to resolve. 
Nevertheless, we estimate for I Zw 18 Af gas ~ 1 x 10 7 Mq, 
on the same order of the 2.7 x 10 7 M© found by van Zee 
et al. <|1998l . The dust mass of IZwl8 is in the range of 
2-5 x 10 3 Mq (Cannon et al. l2U02[) . Therefore, we obtain a 
dust-to-gas ratio ~ 10 -4 , which is roughly consistent with 
the model calculation at the age of ~ 20 Myr. 

4.2.2. SBS 0335-052 

The gas number density is roughly n ~ 10 3 cm~ 3 (Izotov 
et al. 119991 Hunt et al. l20U4a|l . and the radius is r SF — 40 
pc. Eq. indicates that M gas ~ 10 7 Mq. The free-fall 
time is estimated to be iff ~ 3 Myr (Eq. [IJ) . In order to be 
consistent with our previous paper, HHF02 (ip — 1 Mq), 
we adopt a slightly large star formation efficiency, e = 0.3. 
Such a high star formation efficiency can be correct if a 
region is affected by the external pressure (Elmegreen & 
Efremov lTMTI . 

In Fig. [7] (solid line), we show the evolution of gas 
temperature, ionisation degree, and molecular fraction for 
SBS 0335—052 with solid lines. The gas finally cools on a 



timescale comparable to the observationally derived young 
age (~ 5 Myr): the temperature drops to 300 K, the 
medium is kept neutral, and the molecular fraction in- 
creases. As mentioned in Section 12.21 the calculation is 
stopped at t = i gas /2. However, since the gas state rapidly 
changes around t = i gas /2 = 5 Myr, it is interesting to 
show the evolution after t = i gas /2, although neglecting 
the gas consumption becomes a bad approximation. The 
dotted lines continuing after the solid lines show the evo- 
lutions for i > i gas /2. 

The evolutions of oxygen abundance and dust-to-gas 
ratio is also shown in Fig. |S] (solid line). The dashed lines 
shows the evolution after t gas /2 (same as Fig.[7J|. We see 
that the metallicity reaches 1/50 solar (roughly the ob- 
servational value 1/41 solar) around t ~ 6 Myr, which is 
marginally consistent with the age constraint by Vanzi et 
al. (|2000) (<; 5 Myr). We should note that this metallicity 
level is reached after i gas /2, which implies that a signif- 
icant fraction of gas in the star-forming region has been 
consumed. The run-away collapse of cooled gas could lead 
to such a rapid gas consumption. The dust-to-gas ratio is 
- 10~ 4 at t ~ 6 Myr. 

We have adopted a gas mass M gas ~ 10 7 Mq for 
the star-forming region, much smaller than the observa- 
tional estimate ~ 10 9 M© by Pustilnik et al. (|2001j> for the 
whole galaxy. However, because of the lack of resolution 
for H I observation, it is difficult to put a constraint on 
the gas mass in the star-forming region. There is also a 
large uncertainty in dust mass. Dale et al. ij^ QOlfl estimate 
Mdust to be 2400 Mq but they neglect the stochastically 
heated dust grains (in this sense, the dust mass is un- 
derestimated), while Hunt et al. (2001) derive 10 5 Mq, 
but they argue that this is an upper limit because of 
the assumed spatial filling factor of unity. From a radia- 
tive transfer model, which however neglects the effects of 
small grains, Plante & Sauvage ©02) find a dust mass 
of 1.5 x 10 5 Mq. Takeuchi et al. lj2UU5|> adopt the evolu- 
tion model of HHF02, and they explain the infrared spec- 
trum of SBS 0335-052 with M dust - 5 x 10 3 Mq, which 
is smaller than the upper limit by Hunt et al. I|2001J) but 
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Fig. 8. Time evolution of a) metallicity and b) dust-to-gas ratio for the model of SBS 0335—052 and I Zw 18 (solid 
and dotted lines, respectively). The dashed line continuing from the solid line shows the evolution after t gas . The 
dash-dotted line represents the metallicity level of 1/50 solar. 



is larger than the possibly underestimated value by Dale 
et al. <HHJ). 

By using M dust ~ 2x 10 3 -10 5 M and M gas - 10 7 M , 
we obtain a range of the dust-to-gas ratio in SBS 0335—052 
of 2 x 10 _4 -10~ 2 , which is larger than that expected from 
our model. Because of the uncertainties in the dust mass 
and the gas mass, further observations are needed to bet- 
ter constrain these numbers for SBS 0335—052. Recently, 
Inoue ( 2003) also has discussed the dust-to-gas ratio of 
SBS 0335—052 with his model but he has reproduced the 
observational dust-to-gas ratio with t <; 10 7 yr. Nozawa 
et al. (|2003p have proposed a SN dust production larger 
than Todini & Ferrara ; 2001 ;. and their model could solve 
the discrepancy between theory and observation. 

4.3. Possible intermediate active/passive BCDs 

The active and passive modes of star formation are clearly 
two extremes on a continuum. In this paper, we have sin- 
gled out two parameters that determine a galaxy's posi- 
tion on this continuum, namely size r$F an d density n. 
Dust properties distinguish the two extremes, but they 
also must be viewed as part of a continuum of possibili- 
ties. 

Figure 1 of HHTIV indicates that active and pas- 
sive BCDs are clearly separated into two sequences in 
the surface brightness - density relation. However, there 
are BCDs which could be considered as intermediate 
between the active and passive extrema. For example, 
Mrk71 (NGC2366), whose metallicity is about 1/10 Z Q 
(Peimbert et al. ll986[) . has a compact star- forming region 
but a low surface brightness (Hunt et al., in preparation). 
Moreover, no SSCs are found in this galaxy by Billett 
et al. (2002). Another possible intermediate galaxy is Tol 
1214-277, whose metallicity is about 1/25 Z Q (Izotov et al. 



2001). It has a dense (n ~ 210 cm" 3 ; Fricke et al. EUOTJl . 
but rather diffuse star-forming region (size ~ 200 pc; Hunt 
et al. , in preparation) . However, since this galaxy is distant 
(104 Mpc), it is difficult to constrain the size star- forming 
region of Tol 1214-277 even with the spatial resolution of 
Hubble Space Telescope. 

5. Implications for high redshift 

Large numbers of primeval galaxies are expected to ex- 
ist in the high-z universe (z J> 5; e.g., Scannapieco et al. 
2003). Because such galaxies must necessarily be chemi- 
cally unevolved, it is useful to extend our concept of the 
active-passive dichotomy to high z. The "active" mode of 
star formation is particularly relevant for high-z galax- 
ies, because of high gas density <; 1000 cm~ 3 (Norman & 
Spaans 1997) and warm dust content. 

Galaxies become IR luminous on much shorter 
timescales than the cosmic age at z ~ 5 (HF02), even at 
metallicities is as low as I Zw 18 or SBS 0335—052 (see also 
Morgan & Edmunds 2003). We have shown that during an 
"active" star formation episode with no prior chemical en- 
richment, dust produced by SNe II efficiently shields the 
UV radiation. This means that the heating by UV photons 
in the ISM is not sufficient to halt gas cooling. Therefore, 
it is important to trace the high-z star formation activ- 
ity in the IR (or sub-mm) (e.g., Chapman et al. [2001). 
Indeed, large extinction has been shown by some samples 
of high-z galaxies (e.g., Meurer et al. 1999. Stanway et al. 

USUI . 

The "active" mode is also characterised by high dust 
temperature (<; 40 K) because of the high UV interstellar 
radiation field. Totani & Takeuchi (2002) have shown that 
the existence of high temperature dust at z ~ 3 is favoured 
to explain the far-IR cosmic background radiation. Thus, 
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it is important to investigate a potential increase of the 
"active" star formation mode at high z (especially the 
increase of a population with high dust temperature), in 
order to quantify the contribution of such galaxies to the 
cosmic IR background. 

Indeed, virtually all high-redshift star formation may 
occur in the "active" regime. It is not clear whether there 
are high z objects corresponding to the "passive" mode 
of star formation. H2 has been difficult to detect from 
damped Lya clouds (DLAs; e.g., Petitjean et al. [2002 
Ledoux et al. I2003[l . In general, stringent upper limits 
of 10~ 6 have been placed on the molecular fractions 
for DLAs. In passive star-forming regions, the molecular 
fraction is kept as low as ~ 10~ 5 on timescales of several 
tens of Myr. Therefore, some of the DLAs could be "pas- 
sive" star-forming objects. However, the low abundance of 
molecules can also be explained by photodissociation by 
the cosmic UV background radiation (e.g., Hirashita et al. 
2001. 



6. Conclusions 

In this paper, we have theoretically modeled the "active" 
and "passive" star formation modes observed in metal- 
poor BCDs. The "active" mode is characterised by dense 
and compact (n <; 500 cm~ 3 and r$F & 50 pc) star- 
forming regions. On the contrary, the "passive" mode 
takes place in diffuse (n < 100 cm -3 and rgF 100 pc) 
star-forming regions. 

In the dense regions observed in "active" BCDs, the 
gas free-fall time is typically shorter than ~ 5 Myr. Such 
a short free-fall time can enhance star formation activity 
and generate an efficient supply of dust from SNe II. An 
accumulation of dust in such a compact region leads to 
a large dust optical depth, and the region consequently 
becomes luminous in the IR. Even in the environment of 
active star formation, the gas retains the physical condi- 
tions favourable for the gas collapse: cool 1000 K) and 
highly molecular (/h 2 ^ 10~ 2 ). The above characteris- 
tics explain the properties of star-forming regions in the 
BCDs categorised as "active" in HHTIV, especially SBS 
0335—052 (i.e., IR luminous, containing hot dust, rich in 
molecules, compact, and dense). 

In "passive" BCDs, the dynamical time of diffuse re- 
gions is longer than 10 7 yr, and the star formation rate can 
be as low as ^ 0.1 Mq. The increase of dust optical depth 
is milder and therefore "passive" star-forming regions be- 
come IR luminous much later in their evolution, > 10 7 yr. 
Such a passive region has a low ^ 10 -4 molecular content, 
and it would be difficult to detect H 2 in passive BCDs. 

However, the physical state of gas is strongly affected 
by the size and spatial distribution of grains. With a fixed 
total dust mass, if the grain radius is large, the optical 
depth of dust against the UV light is small. This means 
for large grains (~0.01 /im), the gas temperature and ion- 
isation degree tend to be high, and the molecular fraction 
tends to be low. As for the spatial distribution of grains, 
efficient shielding of UV light takes place in a screen ge- 



ometry, which leads to efficient cooling and molecule for- 
mation. On the contrary, in the mixed geometry, because 
of inefficient shielding, the gas is heated to nearly 10 4 K, 
and the molecules are efficiently dissociated. 

We also discussed the evolution of dust-to-gas ratio 
and metallicity. In particular, our model is consistent with 
the observations of I Zw 18 and SBS 0335—052, but future 
observations of dust in metal-poor ([O/H] ^ —1) BCDs 
are needed for further constraints. The consistency also 
implies that around 10-20 % of metals supplied from SNe 
are in the dust phase. 

We have shown how differences in two physical pa- 
rameters of a star-forming region, its size and density, can 
lead to substantially different evolution over time. The dis- 
tinction has been made between "active" and "passive" 
modes, but such a dichotomy is perhaps misleading, since 
they are rather extremes on a continuum. Pressure must 
also drive evolution of a star-forming region, and may de- 
termine a region's initial size and density. However, after 
the onset of star formation, dust shielding of UV photons 
determines the fate of star-forming regions, which finally 
bifurcates into active and passive regimes. 

Acknowledgements. We thank the anonymous referee for use- 
ful comments which improved this paper considerably. We are 
grateful to T. X. Thuan, A. Ferrara, M. A. R. Kobayashi, and 
H. Shibai for stimulating discussions on galaxy evolution. H. 
H. was supported by the Research Fellowship of the Japan 
Society for the Promotion of Science for Young Scientists. We 
fully utilised the NASA's Astrophysics Data System Abstract 
Service (ADS). 



Appendix A: Equilibrium dust temperature 

We derive and discuss the equilibrium dust temperature 
in Eq. I|21|) . We start from the following equation between 
absorbed and emitted energy of a grain: 



CU\jyTTa' 



poo 

: Q UV = Ana 2 / Q m {a, X)irB x (T dust ) dX , (A.l) 
Jo 



where u\jy is the radiative energy in UV defined in Eq. , 
a is the grain radius, i?A(7d U st) is the Planck function 
estimated at the wavelength of A and the dust temperature 
of Tdust, and Qvv and Qm(a, A) are the absorption cross 
sections of a grain normalized by the geometrical cross 
section in UV and IR, respectively. We assume that Quv 
is independent of A (Draine & Lee I1984"|) . and adopt the 
following form for QiR,(a, A): 



Qm{a, A) = 



InAa 



(A.2) 



In the main text, we assume that {3 = 2 (Drapatz & Michel 
H577I Shibai et al. 113991 Takeuchi et al. l2UU3|l . but we 
generalise the wavelength dependence. Since the Planck 
function is written as 



2hc 2 /X 5 



exp(hc/XkT) - 1 



(A.3) 
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where k is the Boltzmann constant, we obtain the follow- 
ing analytic solution for Td us t: 



^dust — 



he 



u\jvQ 



UV 



l/(4+/3) 



k V l6n 2 AahcI(3 + (3) , 
where the function /(a) is dchncd as 



1(a) 



dx (a > 0). 



(A.4) 



(A.5) 



Since 1(5) = 8tt 6 /63, we obtain Eq. for (3 = 2. 
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